Pore-polymer interaction reveals non-universality in forced polymer translocation 
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We present a numerical study of forced polymer translocation by using two separate pore models. Both of 
them have been extensively used in previous forced translocation studies. We show that variations in the pore 
model affect the forced translocation characteristics significantly in the biologically relevant pore force, i.e. 
driving force, range. Details of the model are shown to change even the obtained scaling relations, which is 
a strong indication of strongly out-of-equilibrium dynamics in the computational studies which have not yet 
succeeded in addressing the characteristics of the forced translocation for biopolymers at realistic length scale. 



I. INTRODUCTION 

Polymer translocation has been under intensive research for 
the past decade due to its relevance for e.g. ultra-fast DNA 
sequencing JD-Ht] and many biological processes @]. Most 
of the computational research on forced translocation, where 
the polymer is driven through a nanoscale pore by a poten- 
tial, has dealt with fairly weak pore potential or force and 
assumed close-to-equilibrium dynamics. Recently, hydrody- 
namics has been shown to significantly affect forced translo- 
cation for experimentally and biologically relevant force mag- 
nitudes 0, |8[] . The pore geometry has been noted to have ef- 
fect on the experimental translocation process (see e.g. iH). 
In addition, the effect of the pore model has been discussed 
in iHQii. However, the significance of the pore model in the 
forced polymer translocation has not yet been determined. 

We have recently made a comparison of the unforced and 
forced translocation ifTTIl . In the first case, the process was 
seen to be close to equilibrium, as expected, and accordingly 
its dynamics turned out to be robust against variations in the 
computational model. However, in the second case the ob- 
served out-of-equilibrium characteristics for biologically and 
experimentally relevant pore force magnitudes, made the dy- 
namics more sensitive to differences in the computational 
models (fl ■ Accordingly, we expect the pore model to 
play a significant role in the forced translocation case ll ill . 

So far, we have consistently used a cylindrical pore imple- 
mented by a potential symmetrical around the pore axis d, 
[U[l2l. However, this differs from the typical pore imple- 
mentation where the pore is formed by removing one particle 
from a wall consisting a monolayer of immobile, point-like 
particles |fl3lfl7ll . 

A notable difference is that the latter pore implementa- 
tion results in a driving potential that is inhomogeneous in 
the direction of the pore axis. In the case of lattice Boltz- 
mann simulations these bead-pore models include both the 
square ifTrj [l7tl . and cylindrically lfl8tl shaped pores. Typi- 
cally, the pore-polymer interaction is solely repulsive but the 
effect of attractive interaction has been investigated in a 2D 
model 0. 

In order to investigate how susceptible translocation dy- 
namics is to variations in the pore model we have imple- 
mented in our computational translocation model also the 
bead pore used in IU4I1 . We compare translocation processes 



for the bead pore and our original cylindrical pore. The model 
dependency is expected to be strongest in the pore region un- 
der strong forcing. We also try to evaluate the effect of the 
pore model for moderate forcing. 

When the force bias inside the pore is small, the velocities 
involved are small enough for the hydrodynamics to be ne- 
glected |0| [Till , and the Langevin dynamics (LD) is a valid 
choice for the computational model. Assuming that the poly- 
mer remains close to equilibrium throughout the translocation, 
the control parameter may be defined as the ratio of the total 
external force and the solvent friction, ( = f to t /£ as done by 
Luo et al. in 0- Then the translocation dynamics would 
depend only on (. They reported this to be the case even for 
a long polymer chain driven through the pore with a signifi- 
cant value of total pore force f to t- The scaling of the average 
translocation time with respect to the chain length, r ~ N@ 
would then be determined solely by (. This is in clear contra- 
diction with our earlier observation that the scaling exponent 
/? increases with the force due to crowding of the monomers 
on the trans side |H, [T^l ■ In addition, the effect of the pore 
model along with other model-dependent factors would in this 
case be negligible ifTTIl . 

The increase of j3 with force cannot hold asymptotically as 
pointed out in fTill . since ultimately very long polymers would 
translocate faster with a smaller pore force. From Fig. 1 b) 
in J^] it can be evaluated that this unphysical condition would 
be seen for N > 10 6 for the pore force magnitudes 1 and 10. 
This is comparable to actual DNA lengths of 10 6 . . . 10 9 that 
are far beyond polymer lengths of N < 1000 used in sim- 
ulations. The translocation dynamics then would have to be 
different for A*" > 10 6 , which is in keeping with our earlier 
qualitative description of the forced translocation [HHl]. The 
polymers were seen to be driven out of equilibrium throughout 
the translocation already for modest pore force magnitudes. 
The continually increasing drag force on the cis and crowd- 
ing on the trans side were seen to change the force-balance 
condition, which is directly reflected on the translocation dy- 
namics. This force balance is sure to be different for very long 
polymers. 

The suggested control parameter £, if applicable, would im- 
ply that the forced translocation could be completely char- 
acterized through varying it, which again is in stark con- 
trast to our observation that the forced translocation is a non- 
universal, out-of-equilibrium process for biologically relevant 
magnitudes of pore force HlH. 



Therefore, it is essential to 
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FIG. 1 : (Color online) Schematic depiction of the two pore models, 
(a) The pores are viewed from the trans side along the z-axis. The 
small red circle depicts the cylindrical pore of diameter 1.2a. The 
bead pore is defined by the eight beads each at distance 1.5a from 
the z-axis. The pore beads are drawn with circles using the LJ poten- 
tial cutoff length 2 1//fi as their radius. The light blue area in the center 
of the pore indicates the region where polymer beads have no inter- 
action with the pore beads. In contrast, the cylinder pore model has a 
damped-spring-like potential that acts on particles everywhere inside 
the pore, (b) Side view. The polymer about to translocate (s = 1) 
is drawn as connected dots. The potentials of the two pore models 
differ in both the xy-plane and along the z-axis. 



determine if this contradiction arises from using different pore 
models. 

This paper is organized such that we first describe the poly- 
mer model in Section HT1 and the translocation models in Sec- 
tion |III] The relation between computational and physical pa- 
rameters is discussed in Section [TV] The results are presented 
and discussed in Section [V] We conclude the paper by sum- 
marizing our findings in Section[VTJ 



II. POLYMER MODEL 

The standard bead-spring chain is used as a coarse-grained 
polymer model with the Langevin dynamics. In this model 
adjacent monomers are connected with an-harmonic springs, 
described by the finitely extensible nonlinear elastic (FENE) 
potential, 



U 
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(1) 



where r is the length of an effective bond and i? = 1.5 is the 
maximum bond length. The Lennard-Jones (LJ) potential 
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U LJ = 0, r > 2 1 /V, 



(2) 



is used between all beads of distance r apart. The parameter 
values were chosen as e — 1.2, a = 1.0 and K = 60/er 2 . The 
used LJ potential mimics good solvent. 



III. TRANSLOCATION MODELS 

To perform simulations of polymer translocation in 3D 
we use our translocation model based on Langevin dynam- 
ics Hill . Hence, the time derivative of the momentum of bead 
i reads as 



(3) 



where £, Pi(t), r)i(t), and /(rj) are the friction constant, the 
momentum, random force of the bead i, and the external driv- 
ing force, respectively, f{ri) is constant but exerted only 
inside the pore. For unforced translocation /(r,) = 0. £ 
and T)i(t) are related by the fluctuation-dissipation theorem. 
The dynamics was implemented by using velocity Verlet al- 
gorithm (U. 

In our model the wall containing the pore was implemented 
as a surface on which no-slip boundary conditions are ap- 
plied to the polymer beads. The two different computational 
pore models to be compared can be described as follows (see 
Fig.D: 

Cylindrical pore. The pore, aligned with the z-axis, is of 
diameter 1.2<r(= 1.26) and length I = nb, where b = 1 is the 
Kuhn length of the model polymer, and n is either 1 or 3. The 
pore is implemented by a cylindrically symmetric damped 
harmonic potential that pulls the beads toward the pore axis 
passing through the middle of the pore in the z-direction. The 
pore is nb long, so the total pore force is taken as f to t = Tif, 
where / is the external driving force applied on each polymer 
bead inside the pore. 

Bead pore. An octagonal composition of eight immobile 
particles was used as the bead pore model, as depicted in 
Fig. [TJ This pore model includes a region where a poly- 
mer bead can reside without interacting with the pore beads, 
see Fig. [TJ The effective pore length is approximately b, so 
ftot = /■ Unlike in our model, in lfl4ll also the wall was made 
of immobile particles. We verified that this has no effect on 
the translocation dynamics by reproducing essentially identi- 
cal results on t for the force magnitudes that were reported 
in H. 

In our previous Langevin dynamics simulations with the 
cylindrical pore fill [Till , parameter values £ = 0.73, m = 16, 
and kT = 1 were used for the friction constant, the polymer 
bead mass, and the temperature in reduced units, respectively. 
Hence, for long times the one-particle self-diffusion constant 
was obtained from Einstein's relation as Do = kT/^m sa 
0.086. Time steps of 0.01 and 0.03 were (previously) used in 
the forced and unforced simulations, respectively. In this pa- 
per, to compare the two pore models, we have used m = 1, 
and kT = 1.2 as in [14], unless otherwise noted. Here, the 
time step is typically 0.001. 

The number of beads N is odd for polymers initially placed 
halfway inside the pore and even for polymers having initially 
only the first bead(s) inside the pore, see Fig. Q] Before al- 
lowing translocation a polymer is let to relax for longer than 
its Rouse relaxation time. We register events when a seg- 
ment s — so in the pore is replaced by the segment so — 1 
or sq + 1- The polymer is considered translocated, when it 
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has completely exited the pore to the trans side. Exit to either 
trans or cis side is regarded as an escape from the pore. 



IV. RELATING THE COMPUTATIONAL AND THE 
PHYSICAL PORE FORCE 

The external driving force, called the total pore force, 
hot = nf, where n is the number of polymer beads inside 
the pore and / the z-directional driving force excerted on each 
bead inside the pore, will be taken as the effective pore force 
in our model. The thermal energy is set to kT = 1.2 in re- 
duced units and the length scale is set by the polymer bond 
length b. The magnitude of the effective pore force is then 
determined with respect to thermal fluctuations by comparing 
ftotb with kT. It is noteworthy, however, that since the pore 
length is of the order of the polymer bond length, there exists 
an inevitable computational artefact due to the discrete poly- 
mer model. The external force / is excerted on beads residing 
in the pore, whose number n does not stay constant during 
translocation. In addition this number is different for differ- 
ent pore lengths, so the pore force does not increase strictly 
linearly with the pore length. These effects are, however, mit- 
igated by the fact that n averaged over translocation time is 
constant. 

Correspondence between computational and physical 
length scale can be established by taking the polymer bond 
length b as the Kuhn length for the physical polymer. In SI 
units the bond length for our FJC model polymer can be ob- 
tained as b — 2X p , where A p is the persistence length, 40 A 
for a single-stranded (ss) DNA ll2lll . The total pore force in 
SI units, ftot, is then obtained from the dimensionless total 
pore force, /, by relating ftotb/ ksT = f tot b/kT. Here, ks 
is the Boltzmann constant, and the physical temperature T is 
taken to be 300 K. Thus the effective pore force of f to t = 1 
corresponds to f tot = 1.2 pN for a ssDNA. To relate to experi- 
ments, in the a-HL pore a typical pore potential of « 120 mV 
would correspond to f to t ~ 5 pN when Manning conden- 
sation leading to drastic charge reduction is taken into ac- 
count 122112311 . The used computational pore force magnitudes 
are thus in the physically relevant range JH]. 



V. RESULTS 
A. Forward transfer probability 

Previously ifTTll . we have shown that during the unforced 
translocation the polymer remains close to equilibrium so that 
the forward transition probabilities 

iM S )oc(l-i + ^-p (4) 

derived from the free energy apply in 3D, where 7 = 0.69. 
Including the pore force in the calculation only introduces a 
prefactor to the above Pf(s) not changing its form, i.e. de- 
pendence on s. For this to hold, the translocating polymer 




FIG. 2: (Color online) Forward transfer probabilities Pf as functions 
of the reaction coordinate normalized with polymer length, s/N. 
The data is given for the bead (bd) and cylindrical (cyl) pore. Here 
£ = 0.7 and N is 255 or 256, depending on the polymer's initial 
position. The pore force f to t has the following values from top to 
bottom: /1 = 5.0 (cyl,bd), f 2 = 1.17 (bd), f s = 0.5 (bd). At the 
bottom are Pf for ftot = 0.1 for the bead (bd, distinct □ from the 
solid curve) and the cylindrical (cyl) pore (red) obtained from sim- 
ulations together with the black solid curve calculated from Eq. © 
for the unforced case. For the ftot values 0.1, 0.5, and 1.17 (f 2 up- 
per curve) the polymer was initially placed halfway through the pore 
s = (N - l)/2. For the f to t values 5.0 (both) and 1.17 (/ 2 lower 
curve), the polymer started from the cis side s = 1. The shape of the 
probability curve depends of the pore model (/1), polymer's initial 
position (/2X and changes with the force. 



has to remain close to equilbrium. Indeed, except for a small 
increase due to the force term, the shapes of the measured 
Pf(s > N/2) for the unforced and forced translocation with 
the total pore force f to t — 0. 1 are seen to be similar, cf. Fig. [2] 
Pf for the two pore models do not differ appreciably for this 
small pore force. Hence we conclude that the transition prob- 
abilities can be obtained from equilibrium framework for pore 
force magnitudes of the order of 0.1. 

For the bead pore already for the force f to t = 0.5, the Pf (s) 
curve differs from the equilibrium shape, as shown in Fig. [2] 
Hence, the effect of the pore force can no more be taken as a 
small perturbation to equilibrium dynamics. In keeping with 
this, the transition probabilities P/(s > N/2) for f tot = 1.17 
are seen to depend on the polymer's initial position, here ei- 
ther on the cis side or halfway through the pore. In other 
words, the transition probability from the current state s to 
the next depends on the path through which this state was ar- 
rived at. The polymer cannot then have relaxed to equilibrium 
between previous transitions. 

Using the cylindrical pore, the probability Pf is seen to be 
almost constant and close to one for all s with f to t = 5, see 
Fig. [2] This is in accordance with our previous results that the 
polymer was seen to almost always translocate with ftot — 
3, which means that also in this case Pf(s > 0) would be 
close to one |8y. However, Pf(s) for the bead pore deviates 
significantly from Pf(s) for the cylinder pore when f to t = 5, 




FIG. 3: (Color online) Results from 3D Langevin dynamics with a bead pore model. Keeping ( = ftot/S, fixed, results for three pairs 
of parameter values (/tot,£): (1.17,0.7), (5,3), and (10,6), are shown, ftot = fx, /2, and /3, respectively, (a) Scaling of the average 
translocation time r with respect to the chain length N, r ~ iV 3 . The scaling exponent /3 has values 1.32 ± 0.03 (dotted magenta line), 
1.48 ± 0.03 (solid black line) and 1.52 ± 0.03 (dashed blue line) for the parameter pairs, respectively. The exponents are from fits for 
N £ [32, 64, 128, 256]. (b) Average transition times t(s) = t(s - 1 ->• s) for N = 128. Plots from bottom to top are for f tot = fx, fi, and 
/3, respectively. The transition time profile changes significantly when increasing the pore force eventhough the proposed contol parameter ( 
is kept fixed. 



see Fig. |2 Thus it is evident that the pore force magnitudes in 
the two model pores do not have an exact correspondence for 
large pore force magnitudes. 

The pore force at which translocation is seen to be 
a strongly out-of-equilibrium process is surprisingly low 
for both pore models. It may be compared with the 
random force in the one-particle Langevin equation sat- 
isfying fluctuation-dissipation theorem in 3D, for which 
(f{t)f(t')) = 6kT£5(t - t'). With the parameter values 
£ = 0.7 and kT = 1.2 used in our simulations, this yields 
/~2.2. 

Already at fairly moderate pore force values, a local force 
balance governs the translocation dynamics. This was first 
proposed by Sakaue 12411 . and demonstrated in simulations by 
us ll8i [l2Tl . Considering a concept of 'mobile beads' near the 
pore, we obtained and described qualitatively the scaling rela- 
tion r ~ N l+V ~ a / ftot, where the parameter a taking into ac- 
count the varying number of mobile beads diminishes toward 
zero with increasing pore force |H,[l2l. Sakaue has recently 
given a fairly quantitative out-of-equilibrium formulation for 
the scaling r ~ N 1+lJ / ftot lHHl- However, neither description 
takes into account the crowding of monomers on the trans side 
shown to be significant for these pore force magnitudes [8|]. 

Here, we present results from simulations of the polymer 
translocation covering the relevant range of the pore force. 
When ftotb/kT S; 0.1, Pf(s) attains the equilibrium shape 
and when ftotb/kT <J 5, P/(s) is close to one. 



B. The proposed control parameter 

Using the bead pore, we first measured the scaling exponent 
f3 for different values of pore force f tot while keeping the pro- 



posed control parameter ( = f tot / £ fixed. For the parameter 
pair (ftot, values (1.17, 0.7), (5, 3), and (10, 6) were used, 
see Fig. [3](a). Similar average translocation times r were ob- 
tained for large pore force magnitudes, f to t — {5, 10}. How- 
ever, t differs appreciably for f to t = 1-17, resulting in a dif- 
ferent value for the scaling exponent (3. 

The measured average waiting (or state transition) times 
t(s) = t(s — 1 — y s) are compared in Fig. [3] (b). The 
shown three plots for t(s) are for the above-mentioned val- 
ues of (ftot,0- I n tne limit ftotb/kT <§C 1, the average state 
transition times t(s) were seen to be similar. This is in accord 
with the transition times in the unforced translocation obeying 
a close-to-Poissonian distribution |0. When ftotb/kT > 1 
the friction of the solvent dominates over the stochastic term 
in the Langevin equation, Eq. ([3j. The transition time profiles 
for the bead pore are qualitatively similar with those obtained 
for the cylindrical pore with different parameters lfl2tl . 

From the different scaling of the translocation times ob- 
tained for large and small pore force magnitudes and the dif- 
ferent average state transition times we conclude that C which 
was kept fixed cannot be a universal control parameter for 
forced translocation regardless of the pore model. This is 
in keeping with our previous finding that the forced translo- 
cation is a highly out-of-equilibriun, non-universal process 
for the biologically and experimentally relevant pore force 
range §,[11. 



C. Pore model matters in forced translocation 

Originally, the need to compare the cylindrical and bead 
pore models arose from the differences in the scaling of the 
translocation time r with the polymer length N in forced 
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FIG. 4: (Color online) Scaling of the average translocation time with respect to the chain length, r ~ TV 3 . The results are from 3D Langevin 
dynamics simulations with both cylindrical (cyl) and bead (bd) pore models with different pore lengths / = nb. (a) The constant total pore force 
ftot = 5. The exponents are from fits for TV G [32, 64, 128, 256]. For £ = 0.7, the scaling exponent P = 1.26 ± 0.02(D), 1.27 ± 0.02(a) 
and 1.36 ± 0.03(O), for Z = 3 (cyl), I = 1 (cyl) and I = 1 (bd), respectively. For £ = 3, the scaling exponent /3 = 1.39 ± 0.02(O) 
and 1.48 ± 0.02(D), for I = 3 (cyl) and / = 1 (bd), respectively. Increasing the bead mass from 1 to 16 increases the absolute value of r, 
and P = 1.40 ± 0.03(«) for £ = 0.7 (cyl). See text for details, (b) For f to t = 0.5, £ = 0.7 (T: bd) the polymers are initially halfway 
through the pore. We obtain P = 1.25 ± 0.04 (green dashed line) with TV G {31, 63, 127, 255, 511}. For ftot = 0.1, £ = 0.7 (D: cyl, 
A: bd) polymers are also initially halfway through the pore. The absolute values of r differ slightly for the two the pore models with shorter 
polymer chains but not so for TV ~ 511 when all chains escape to the trans side. For a very low force ftot = 0.01(O : bd) even the longest 
(TV = 511) chains may escape to the cis side. We obtain /? = 2.2 ±0.1 (solid black line). For reference, results with ftot = 0.1, £ = 6 (o: 
bd) are shown. Inset: The scaling exponent P as a function of the friction constant £, f to t = 5 (bd). For £ = {0.7, 1.5, 3.0, 6.0}, we have 
/3 = 1.36 ± 0.04, 1.44 ± 0.03, 1.48 ± 0.03, and 1.50 ± 0.03, respectively. 



translocation obtained using these two pore models. The 
translocation times averaged over at least 1000 and at most 
2500 runs are shown in Fig.|4](a) as functions of TV. The total 
pore force is constant, ftot = 5, and chosen to be within the 
experimentally relevant range (8[ ■ For the friction parameter £ 
three values, 0.7, 3, and 6, are used. For the two separate pore 
models we find clear differences in both the absolute value of 
r and its apparent scaling. Regardless of the value for £, we 
obtain a smaller scaling exponent B for the cylindrical than for 
the bead pore. In addition, increasing £ from 0.7 to 6 increases 
B for both pore models, see inset in Fig. |4](b). This in part 
explains the different results obtained with the two pores as 
the effective friction experienced by a chain is different inside 
them. This difference is naturally enhanced at higher veloci- 
ties. 

Our previous results were obtained with the polymer bead 
mass m = 16 due to the requirement of the stochastic rotation 
dynamics that the polymer beads should be heavier than sol- 
vent beads JH[T2l. Therefore, we checked how the bead mass 
affects the characteristics of the translocation time. Increas- 
ing the bead mass from 1 to 16 is seen not only to increase 
the average translocation time r but also to change its appar- 
ent scaling with polymer length r ~ TV 3 , see Fig. [4] (a). To 
be certain, we checked this using two independent Langevin 
algorithms. 

In accordance with our previous findings for cylindrical 
pore, j3 was seen to increase with pore force also for the 
bead pore. For f to t = {0.5, 1.17, 5, 10}, the scaling exponent 
(3 = 1.25 ± 0.04 (Fig.E(b)), (3 = 1.32 ± 0.03 (Fig. 0(a)), 



8 = 1.36±0.03(Fig.]4](a)), and (3 = 1.37±0.03(not shown) 
were obtained, respectively. Surprisingly, for the cylindri- 
cal pore B, which increased substantially with increasing f to t 
when m = 16 0, [lH, did not show such strong tendency 
when m = 1. Changing m seems then to change even the 
qualitative characteristics of forced polymer translocation. 



The change of (3 in the scaling r ~ TV^ due to a change 
in £ is hardly surprising. However, B changing with m is 
a bit more subtle. It suggests that we are actually looking 
at the forced translocation in a continually changing or tran- 
sient stage. This is in accord with our previous finding that for 
polymers of lengths that are used in simulations, TV < 1000, 
the number of moving polymer beads changes throughout the 
translocation and thus continually alters the force balance con- 
dition resulting from the forced translocation taking place out 
of equilibrium iH . Also the crowding of the polymer beads on 
the trans side, whose relaxation toward equilibrium is slower 
than the translocation rate changes the force balance continu- 
ally. Hence, the simulated forced translocation processes do 
not reveal the asymptotic (TV > 10 6 ) characteristics for ex- 
perimentally relevant pore force magnitudes. Then reporting 
scaling exponents for such forced translocation seems unwar- 
ranted. The remaining question of interest then is, why does 
the forced translocation exhibit the scaling albeit with varying 
exponents? 
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D. Translocation at low pore force 

For simulations at low pore force the chains were initially 
placed halfway through the pore unlike at large pore force 
where the polymer was placed so that only its end was inside 
the pore. At the very low force of f tot = 0.1 the pore model 
was found to affect only slightly the absolute value of r, see 
Fig.|4](b). For the bead pore model chains of length N < 127 
had a finite probability to escape to the cis side. Unlike for the 
bead pore, even some of the 255 beads long polymers escaped 
to the cis side for the cylindrical pore. This is probably due 
to the cylindrical pore aligning the polymer chain toward the 
pore axis in the middle thus effectively reducing the friction 
between the polymer and the pore. 

For the pore force f to t = 0.01 we obtain j3 = 2.2 ± 0.1, 
see Fig. [4] (b), which is the expected exponent for unforced 
translocation, (3 = 2v + 1 = 2.2, since v = 0.6 has been 
measured for the swelling exponent in our model lfl2ll . Ap- 
proaching the unforced case makes the translocation dynam- 
ics increasingly robust to variations in the model, which is a 
natural consequence of the small pore force not dominating 
over entropic forces. We regard the pore force magnitude 0.1 
already small enough for the translocating polymer to remain 
close to equilibrium, see Fig. [2] However, we did not obtain 
a clear scaling for f tot = 0.1 for either pore models, but (3 
close to the unforced translocation value 2^ + 1 was obtained 
for chains shorter than N = 127 while a lower value of (3 
was obtained for longer chains. This applies for both £ = 0.7 
and £ = 6, see Fig . El (b). We expect the lower value of j3 to 
be 1 + v (see Il26l0 . although a wider range of N would be 
needed to confirm this. This apparent cross-over behavior for 
low pore force is outside the scope of the present paper and 
calls for a separate, more thorough investigation. 

For the presumably small pore force ftot = 0.5 we ob- 
tained f3 = 1.25 ± 0.04. This differs considerably from 
P = 1.58 ± 0.03 reported in (H. From Fig. |2] it can be 
seen that Pf for f to t = 0.5 clearly deviates from the close- 
to-unforced form of Pf for f to t = 0.1. Together with the 
low value of (3 this observation suggests that f tot = 0.5 is 
large enough to drive the polymer more and more out of equi- 
librium as the translocation proceeds. The different values 
for (3 here and in flill also support this view, since the two 
simulations started from different initial positions. When the 
relaxation time of the polymer toward thermal equilibrium 
is slower than the translocation time, the process possesses 
memory, or correlation over time, i.e. the memory function 
M(t — to) 7^ S(t — to). Thus, a polymer starting half-way 
through the pore is in a different state than a polymer that has 
arrived at this position and started from another initial posi- 
tion. 

The data presented for different pore force magnitudes in 
Sections [VCl and IV Dl show that for f to t > 0.5 the value of j3 
increases with increasing pore force JH [HI] also for the bead 



pore. For very low pore force f3 is found independent of the 
pore force. 

VI. SUMMARY 

In summary, we have simulated forced polymer translo- 
cation in 3D by using Langevin dynamics. We have imple- 
mented two pore models in our algorithm, (i) a pore sur- 
rounded by eight immobile beads, which we call bead pore 
and (ii) a cylindrical pore where a damped harmonic poten- 
tial confines the beads inside the pore region. The present 
study compares the effect these two pores have on the forced 
translocation. 

We measured the forward transition probabilities Pf (s) to 
determine the range of pore force that would be sufficiently 
large to include cases where polymer translocation takes place 
close to and strongly out of equilibrium. We found that the 
polymer remains close to equilibrium when the total pore 
force ftotb/kT ^ 0.1. Here the translocation dynamics was 
found to be robust to variations in the translocation model, 
such as the details of the pore. Pf(s) was found to deviate 
significantly from the close-to-equilibrium form already for 
a pore force as small as 0.5. The polymer was found to be 
driven far from equilibrium when f tot b/kT 2; 5. 

For small pore force magnitudes the forced translocation 
processes are identical for the two pore models. However, 
the translocation characteristics were found to be increasingly 
model dependent when the pore force is increased. This is a 
natural consequence of the dynamics of the forced translo- 
cation being determined by a continually changing force- 
balance condition when the pore force is large enough, i.e. 
{ftotb/kT 2 0.5). Accordingly, it seems that universal ex- 
ponents for the forced translocation cannot be found in the 
biologically relevant pore force regime. In addition, attempts 
to define a control parameter whose magnitude would con- 
sistently determine these exponents would seem futile, which 
we showed for one proposed candidate by using the bead 
pore. Qualitatively the forced translocation exhibited simi- 
lar characteristics with both pore models, most notably the 
increase with the pore force of the exponent determining the 
relation between the average translocation time and the poly- 
mer length. In our view this is another indication of the contin- 
ually changing force-balance condition governing the highly 
non-equilibrium forced translocation process. 
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